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■ Abstract 

' A scheme for the generation of non-classical pairs of photons in atomic 

^ I vapours is proposed. The scheme exploits the fact that the cross correla- 

O^. tion of the emission of photons from the extreme transitions of a four-level 

cascade system shows anti-bunching which has not been reported earlier 
and which is unlike the case of the three level cascade emission which 

^ ■ shows bunching. The Cauchy-Schwarz inequality which is the ratio of 

' cross-correlation to the auto correlation function in this case is estimated 

^Nj , to be 10'^ — lO'' for controllable time delay, and is one to four orders of mag- 

' nitude larger compared to previous experiments. The choice of Doppler 

I free geometry in addition to the fact that at three photon resonance the 

. excitation/deexcitation processes occur in a very narrow frequency band, 

' ensures cleaner signals. 

o 

1 Introduction 

' Non-classical sources of light include triggered single photon emitters [T] as well 

^ I as photon pair emitters [5] where the first photon heralds the arrival of the 

second one. The recent elegant experiments on the generation of non-classical 
photon pairs have successfully reported the temporal control of emission of an 
anti-stokos photon by programmable time delay. Polyakov et al [3], report a 
violation of the classical Cauchy-Schwarz inequality in the photon correlations, 
and they further report that the ratio of cross-correlation to auto-correlation 
j^max _ ± 57. All these experiments have been done on cold atoms in 
magneto-optical traps (MOT) . Recently, this experiment has been repeated in 
atomic vapours at room temperature [4] and the results are comparable with 
the MOT results. While in the MOT set up the magnetic and quadrupole 
fields induces decoherence in the correlations, in the latter case the inefficient 
absorption/excitation due to Doppler broadening is the limiting factor. 
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In this paper I propose a scheme for the generation of non-classical pair of 
photons wherein Doppler effects are minimised even at room temperatures, and 
the degree of violation of the Cauchy-Schwarz inequality is very large ~ 10^—10^. 
The key feature we exploit here hinges on the fact that the cross correlation be- 
tween the emission of the extreme transitions in a four-level cascade system 
shows anti-bunching, which is in contrast to the well known fact that cross- 
correlation of the emission from three level cascade systems shows bunching 
O [6] . In addition, one can obtain very strong correlations between the pho- 
ton emissions from the extreme transitions by controlling the driving fields. To 
achieve this I consider here a model which consists of the 4-level cascade system 
in Rubidium which is an extension of the well studied 3-level cascade system in 
the context of the absorptive as well as the dispersive property in electromag- 
netically induced transparency (EIT) systems [ZlinilH]. In addition to the 3-level 
cascade system consisting of 5si/2,5p3/2 and considered in the previous 

examples, an additional level, say, one of the other hyperfine levels of 6^3/2 
may be included and a strong coupling introduced between different hyperfine 
levels of 5p3/2 through a radio-frequency (rf) field (see Figl). The role of this 
'sandwich' coupling is to change the usual 'bunched' emission from the cascade 
system to antibunched emission, with the difference that the cascade emission 
is interrupted by the rf field after the first emission. This control field helps in 
inducing strong correlations between the uppermost and the lower most tran- 
sitions. It further provides us with a handle to control the time delay between 
the successive emissions of photons from the extreme transitions. The combined 
feature of the violation of two classical inequalities namely (i) antibunching and 
(ii) Cauchy-Schwarz inequality in this scheme makes it very attractive. Fur- 
ther the problem of the absorption/emission into other frequency modes close 
to resonance (due to Doppler broadening) which leads to the weakening of the 
signals [4] is circumvented in the present scheme because of the intrinsic narrow 
absorption features due to atomic coherence effects induced by three-photon 
resonance. 

The plan of the paper is as follows: In section II the model under consid- 
eration is described in detail and the equation of motion for the atom-|-field 
system is set up. In section III the relevant photon correlation functions are 
evaluated numerically and in some special cases analytical expressions for the 
correlation functions are presented, which are evaluated perturbatively. Further, 
these results are compared with the exact numerical results for the correlation 
functions. 

2 The model 

The system considered here consists of an ensemble of four-level ladder atoms 
(Figl). This could correspond, eg., to 5si/2, two hyperfine levels of ?>p^/2 and 
5c?5/2 of Rb TO]. I would like to recall here that this was exactly the scheme 
used by Banacloche et al 7J and Fulton et al [9] for studying EIT. Labeling 
the energy levels 5si/2, 6^3/2 as |1 > and |2 > respectively, level |3 > may 
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Figure 1: Four- level system interacting with three driving fields of strengths 
Qi, flrf and their respective detunings Ai, A^/ and A2. r2,r3,r4 are the 
decay constants of the corresponding levels. 

be identified as one of the hyperfine levels of 6^3/2 (F=0,l,2,3 of ®^Rb and 
F=l,2,3,4 of ^'■'Rb ), with level separations typically of the order of lO'^MHz. 
Level |4 > may be chosen to be 5^5/2- The only dipole allowed transitions are 
|1 ><-^ |2 >, |2 ><-^- |3 > and |3 ><-^- |4 > in the limit of the hyperfine splittings 
being larger compared to the Rabi frequency [18] . Comparing this scheme with 
the 3-level cascade system of references [3 [HI [S], levels |l>,|2>,|3>of these 
setups correspond respectively to levels |1 >, |2 > and |4 > of the present model. 
In addition we have introduced another hyperfine level which we have labeled 
as |3 >. 

For this system, the two transitions |1 >— + 2 > and |3 >— + 4 > respectively 
may be chosen to correspond to wavelengths 780nm and 775. Snm, as in the 
schemes used in references [H El [8] . Apart from these two driving fields in the 
optical region, we consider yet another coupling between the hyperfine levels 
( belonging to 6^3/2) |2 > and |3 > through a rf field of strength Qrf- Such 
couplings between the hyperfine levels have been considered in the study of ab- 
sorption in 4- level systems [T^]. The strengths of the optical fields are denoted 
by f2i and fis, and the respective detunings by Ai and A3. The corresponding 
decay parameters are [TO] r2,3 = 6MHz, r4 « lMHz{0.97MHz). To make the 
three photon interaction Doppler free we choose the two beams in the optical 
region, as already mentioned, to be counter-propagating. They couple respec- 
tively transitions |1 ><-^ |2 > and |3 >^ |4 >. The optical frequencies, being 
nearly equal, will give null contribution to the linear Doppler shift. While, the 
linear Doppler shift of the rf field (in the atomic rest frame) at this frequency is 
very small and hence negligible. This choice for minimizing the Doppler effect is 



preferable compared to the phase-matched geometry used in four-wave mixing 

as this offers a larger interaction volume compared to the latter. 

The master equation describing the interaction of the four-level atomic sys- 
tem with three driving fields in the rotating wave approximation is given by 

dpi2 



dt 

dp23 
dt 

dp34 

dt 

dpi3 
dt 

dpii 

dt 

dp24 
dt 

dp22 
dt 

dp33 
dt 

dpAA 
dt 



(-zAi - r2/2)pi2 - i0.l{p22 - Pll) + i^rfPl3 

{-iA2 - {T2 + T^)/2)p23 - ifiipis - 

i^r{{P33 - P22) + i^3P2i 
{-iAs - {T3 + r4)/2)p34 - inrfP24 

-i^siPM - P33) 

(-i(Ai + A2) - r3/2)/)i3 - ifliP23 
+iflrfPl2 + i^SPli) 

(-i(Ai + A2 + A3) - Ti/2)pu - i^ip24. 
(-i(A2 + A3) - (Fa + Vi)/2)p2i - i^ipu 

-iflrfPSi + i^3P23) 

-T2P22 + i^i{P2i - P12) 

+iflrf{p23 - P32) + 723P33 + 724/344 
-r3P33 + i^3{P3i - Pis) 
-«^r/(P23 - P32) + 734P44 
-T4P44 - iQ3{p34 - P43) 



where 7jfe are the transition rates from the level | A; > to |i >, A, = Wj^i+i — are 
the laser detunings [pij = pji* and Trp =1). Throughout this paper we assume 
7jj_l_i = 1 and the rest to be zero. I obtain analytical solutions by perturbatively 
solving these equations in the Laplace space in the special cases where the rf 
field strength is either i) large, or ii) small compared to the other two driving 
fields. 



Ccise (i) Large flrf 

In this regime we assume that flrf > l^i, ^3. Treating the fields fli, z = 1, 3 per- 
turbatively, and labeling Vi = P12, tp2= P23, V's = P34, V'4 = Pi3, V'5 = Pi4, tpe = 
P24, i>7 = P22,i>8 = P33, i>9 = P44 and '^i'i(s) to be the Laplace transform of tpi{t), 
the equation of motion in the Laplace space assumes the simple form 

(s + f2)4^^(s) = V7(0) - 20^/7mV;f + 
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{s + f 3)4'^(s) = MO) + flm4^\s) - 

{s + f4)4^\s) = V9(0) + 203/mVii'^(s) (1) 

where Slrf is treated upto all orders and the fields Oi,3 are treated perturba- 
tively. For simplicity I have assumed A^/ = A; = 0, i = 1, 3 and replaced Ti/2 
by Ti everywhere. Further, the first order solutions for ip[^\ ip^^ndtpo'^ 
are obtained by solving 

(S + f 2)^5^^ (S) = ifir/^i'^ (s) + ini{24"^ (s) + 

-inrf4'\s) (2) 

The solution for Pii{t) depends on the choice of the initial condition which will 
be discussed in section III. 

Case (ii) Small flrf. 

I next consider the case when the rf field is weak. Considering the contribution of 
f^i,3 upto all orders and treating Qj-f perturbatively,the second order equations 
are given by 

(s + f 2)Vif ^ = 20i7mV;f ' - 2nrflm4^^ + 

Vif +V7(0) 

(s + f3)4^^ = -2n37mVi3(2) + 2nrflm42^ + 

+ ^-8(0) 

(s + r4)4'^ = 2r!3 ^ + V'g (0) 

(S + Fg + f 4) V^^'^ = -i^rf4^^ - (4'^ - ) 

{s + f 2)vif ^ = -ij^i(2v;f ) + + 4^^ - 1) + 

iftrf44^ (3) 

where the first order ■0i are obtained by solving 

(s + f 3 + T2)4^^ = -ini44^ + in34^^ 



5 



-zl7,/(V'8(0)-^7(0)) 

(s + f 2 + f 4)v^^^' = -miij^^^^ + iVL^ii,''^^ 

{s + f 4)4" = -ini^^'i^ + in^^^p (4) 

The solutions of this equation are discussed in the next section. 

3 Correlation Function 

Our interest is in determining the second order correlation function G\j (ri, ii, r2, 
which gives the correlation between the fluorescence signals li and Ij at times 
ti = t — ri/c and t2 =t — r2/c, andri,r2 are the positions of the detectors. Here 
i,j=l,3 corresponds to the emission of the g.s. and the upper excited state re- 
spectively. The correlation function is determined by evaluating the expectation 
values ofthe product of two-time operators < E~ {ri,ti)E~ {r2,t2)Ej' {r2,t2)E^ {ri,ti) > 
for i=l and 3 (since we are interested only in the uppermost and the g.s. tran- 
sitions) where E^ and E^ are the positive and negative frequency components 
of the field emitted. 

The normalized temporal intensity correlation function {ti,t2) can be 
written in terms of the atomic operators in the far-zone approximation [131 114j 
to be 

(2),. <cr+ih)a+it2)aj{t2)a-{h)> 

5y (11,^2)- +_ (5j 

< ala^ >ss< crjaj >ss 

for unit detector efficiency, where < cr^cr^ >ss corresponds to the steady state 
signal of the ith mode, = |« -I- 1 >< i\ and it~ = \i >< i -I- 1| being the 
atomic operators. The two-time expectation values are evaluated using the well 
known Onsager-Lax Quantum regression theorem [17] which states that if for a 
complete set of system operators Af^, /i = 1, 2...n the one time expectation value 
satisfies 

< A^it) J2 ^f^At, t') < A^it') >,t' <t (6) 
then the two-time expectation value takes the form 

< d„{t')A^{t)dp(t') E.M^At,t') 

<do.inA,{t')d[}it')>, 

t' < t (7) 

for any two system operators Oa and Ofs. Essentially the correlation functions 
(two-time averages) satisfy the same equation of motion as that of the expec- 
tation values of one-time averages [M] . The c-number coefficients A^^jiy are 
derived from the solution of the Heisenberg equation. Using the operator alge- 
bra aijUki = Sjkcru and Eqn.(5), the two-time averages for the four-level cascade 
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(2) (2) 

Figure 2: (Color online) The cross-correlation functions a) g^-^'{T), b) 532 (t), 
and c) g^iir) function of 7r with fi^ = fis = 47, and flrf — 2O7. 

system may be expressed as (TSl HH E] 

ffn W = < P22{t) >p(o)=|i><i| /P22 
933 (t) = < pair) >p(o)=|3><3| Ml 

9^31 (t) = < P22{t) >p(0) = |3><3| I P22 (§) 

wherer = ^2—^1- Here pff denote the steady state values. The functions ^11 (t) 
and (t) ^-'^"^ ^^'^ auto-correlation functions of the g.s. emission (|1 >^ |2 >) 
and the upper excited state emission (|4 >^ |3 >) respectively. The cross- 
correlation function 331 (t) gives the conditional probability for the emission 
of a photon from the |1 >^ |2 > transition at a time r given an emission of 
a photon from the |3 |4 > transition at time t = 0. Here we have used 
the notation of reference :14j: < Pii{T) >p(o)=|j><j| denotes the population 
pa at a time r given the initial condition that the atom is prepared such that 
p(0) = \j >< j\. The time dependent functions Pijir) in the general case 
is determined numerically while in some special cases I provide the analytical 
solutions. I first consider the case when flrf is large. Solving the eq.(l) and (2) 
in the Laplace space wherein both fli and fl^ are treated perturbatively while 
flrf is considered upto all orders , the correlation functions are given by 

9?i\r) =l + ELi«.e-"'- 

gg\T) = 1 -I- e-^^^bo + &ie-"i^ -I- 626-"^^) + 

933 (t) = 1 + e-^^^ico + cie-"i^ + cae-"^'') 

(9) 
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Figure 3: (Color online) The time delay between the emission of the photons 
from the extreme transitions as a function of the driving field strengths i) f2i 
with Q2 — 127 and fls — 47, ii) f22 with f2i = f23 = 47, iii) ^3 with fli — fl2 = 
47, 
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where Oj = aija^ h = 6j/6jvandci = Ci/cN is the normaUsation factors 
aN,bN,CN . The coefficients are given by 



UN = 




bN = 


2fi2^ni/f4(nti(aO) 


Cn = 




cq = 


-ao/f4(f4 - Q;i)(f4 - a2) 


Cl = 


— Q!o/a2(ai — oi2){T4 — 0:2) 


C2 = 


ao/ai{ai — a2)(r4 — ai) 


ao = 


20i(f2 + f4) 



The roots in the exponent arc given by the conjugate pairs 



ai,2 = -r2 - f 3 ± i<^; <^ = ^4fi2^ - (fa - rs)^ (10) 
and the cubic roots 

a3= -|j(f2 + f3)-4J]2^. + f 

a4,5= -^(f2+f3) + 12(l±iv^)fi2^-i^ 



where 



a= 33(3^2 (fa + fa) + 



i^324(f 2 + f 3)'^^^/ + 6912^^6^)3 
Note that the coefficients c, satisfy the condition 



^ Ci + cjv = 



i=0 



The coefficients ai,bi,i = 1,5 arc functions of a-i and since the expressions are 
very lengthy we have not hsted them here. However, the solution for pair) 
in the Laplace space which are displayed in the Appendix, yield an important 
identity amongst the coefficients 



ajv + ^ Oi = 

5 

bN + Y,^i= 



i=0 



While the conditions for the coefficients o, and Cj reflect the antibunching nature 
of the auto-correlation function, the condition for bi implies the anti-bunching 
nature of the cross correlation function g^^{T). This contrasting feature which 
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is exhibited by a 4-level ladder system is unlike the behaviour of the cross 

correlation fmictions (t) in a 3-lcvcl cascade system which is non-zero at 
T = . A comparison of the various cross-correlation functions of a 4-level 
cascade system will be discussed in detail shortly. 

The expressions in Eqn.(9) for the correlation functions indicate that the 

(2) (2) ~ ~ 

decoherence of 533 (t),.?.^^ (r) is controlled by the decay constants r2,r3 and 
r4 while the dephasing of the g.s. autocorrelation function is due to f 2 and f 3 
only. The imaginary part in the exponents indicate that the the strength of the 
correlation oscillates with f7,.y . 

In the limit of weak ^1^/ the solutions for the correlation function are given 

by 

g[f'{T) =l-e-f='"(Cos(20ir)- 

diS'm(20ir)) 

g^^{T) =1 + e-^^'^{xo + xiSin{fliT) + X2Cos{fliT) + 

di'^ir) =l + E-=oyie-^^+3r (11) 

where di = f2do/2f2i, with do = 2fif/(4ri^ + f|). The roots ai,i = 3,5 are 
obtained by replacing f2 ^ f3,f3 — > T^&ndflrf i^s in Qj. Further 

a6,7= -(f 2 +f3 + f 4) T V^3 + ^4>i(t>2 
as,9= -(f 2 +f3 + f 4) T - 20102 

where 

03= (4(fi? + fi3)-(fi+fi + fl-2f3f4) 

01= V(40f-ri) 
02= V4f23 - (rs - r4)2 

In this regime, the g.s. autocorrelation function resembles the second order 

correlation function of the usual two-lcvcl system. This is intuitively obvious 
since the upper level coupling is weak. Hence the dephasing of this function 
depends on r2. The strength of the correlation oscillates with Oi as shown by 
equation (10). While the dephasing of 533^ (r) and g^i (t) depends on Fj, i = 2, 4. 

The strength of the correlation now depends on fii and ^3. 

f 2) 

It follows from the Eqn.(lO) that 5ii (t) = for r = 0. Again, it may 

be verified using the solutions for < p^^ >p(o)=|3><3| and < P22 >p(o)=|3><3| 
in the Laplace space given in the Appendix, that the auto-correlation function 

(2) (2) 
(733 (r) and the cross- correlation function g^i (t) are zero at r = . Since 

the expressions for Xi and are too lengthy, we list all the solutions in the 

Laplace space in the Appendix. Thus, in cither of the regimes, I would like to 

(2) 

emphasize that the function 531 (t) = at t = 0. This is further evident from 
the numerical results which are solved exactly and valid for all strengths of the 
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coupling fields (Fig4c). Note that the cross-correlation is very strong for jt < 1. 
Whenever the emission of the two photons occurs for time delay smaller than 
the lifetime of the atom, it necessarily means that both the photons belong to 
the same cascade emission. Hence the photons show very strong correlation for 
times smaller than 7~^. For larger time delay between the emissions of the two 
photons decoherence weakens the correlation. 

Comparison of various cross- correlation functions 

For easy comparison with the 3-level cascade system, we consider the most gen- 
eral 4-level system and label the emission from the >^ \i > as the 
mode. The cross-correlation functions glj {T),i ^ j which are proportional to 
^ Pj+i-j+ii''') >p(o)=|i><i| listed in the Appendix, give the probability of emis- 
sion of photon in the j*'' mode at a time t given the emission of a photon in the 
i*'' mode at t=0. The function g^^ir) corresponds to the usual cross-correlation 

(2) 

function of the 3-level cascade system and is nonzero at r = 0. Again g^2 {''') 
is also nonzero at r = 0. This is easily understood by looking at the solu- 
tions in the Laplace space (see Appendix). Both 532 (t) and g2i {t) which are 
respectively proportional to < pssir) >p(o)=|3><3| and < P22(t) >p(o)=|2><2| 
have contribution from terms like s^/(s + a3)(s 4- a4)(s + a^) in the Laplace 
space the inverse transform of which has the form ^iC^"''^ where Ai = 
af/{ai-aj){ai-ak), i 1^ j 1^ k,i,j, k = 3,5. Clearly, at r = 0, the J2i=3 7^ 

. I would like to recall here similar results due to Loudon [5] for a 3-level cascade 
(2) 

system, where g2i (t) is nonzero at r = . On the other hand, the distinguish- 

(2) 

ing feature 531 (t 0) = is due to the absence of such terms. Thus, while 
(732 (r) and g2i (t) both show bunching, g^i {t) is zero at r = 0. The exact 
numerical solutions which are valid for all strengths of the coupling fields, is 
illustrated in Fig2. 

To understand the underlying mechanism behind this feature consider the 

limit flrf — > and f2i,i = 1,2 ^ 0. In this limit the correlation function 

g^iir) is proportional to e^'^'"^/^)'^. i. e. the correlation for the two photon 

emission pathway |3 >^ |2 >^ |1 > depends on e^^'"^/'^^'^. Likewise , the 

(2) 

cross-correlation for the emission pathway |4 >^ |3 >^ |2 >, 1732 (t) is propor- 
tional to e^''"^/^)'^ and hence both these correlations are nonzero at r = 0. On 
the other hand the gfi^(T) = 2/(r3 - r2)(e-(^^/2)r _ g-(r3/2)T)^ rpj^jg implies 

that the cross correlation of the photons emitted by the extreme transitions 

gets contribution from both the two-photon pathways. At r = this function 

vanishes. In other words, there is an exact cancellation of emission due to the 

two pathways |4 >— > |3 |2 > and |3 >-+ |2 >^ |1 > at r = 0. Thus the 

(2) 

antibunching nature of 331 (t) is intrinsic to the four-level cascade system and 
is independent of the driving field strengths. Since this implies a finite time 
delay for the emission of the second photon and since the offset is provided by 
the strength of the driving fields (Fig4c), in principle one has a handle on the 
control of this time delay r^. Fig3 shows the variation of with = 1,3 
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where Td is the time delay at which the cross-correlation peaks. It is clear from 
Fig3 that larger the field strengths, smaller the time delay between the emis- 
sions. As ^2 or ^3 is increased the uppermost level gets populated much faster 
and hence the peak of the cross-correlation function occurs at earlier Td ( This 
is reflected in Fig4c (ii) and (iii) for ilrf — IO7 and 2O7). While increasing fii 
increases the population in level |2 > till saturation. Since this does not affect 
the cross-correlation significantly, the variation of Td with Qi is not as significant 
as with ^2,3- The control will be better with the use of pulsed excitation with 
programmable time delay (of the order of a few 7's) between successive excita- 
tions [2 [3] . This provides an advantage over the other schemes like the photon 
down conversion or the usual three level cascade emission where the emission of 
both the photons is almost simultaneous (bunched) and hence there is no way 
of delaying the emission of the second photon after the first photon is emitted. 

I would like to recall here that antibunching is the violation of the classical 
inequality: ^^^•'(O) > 5'^'(t). This inequality is satisfied by 'classical fields', i. 
e., the Glauber-Sudarshan phase-space function corresponding to the states of 
the field has positive distribution. Violation of these inequalities would imply 
that the distribution function corresponding to the state is not well behaved 
(not positive definite) and are non-classical. Having demonstrated the anti- 
bunching nature of 531 (r) in this section, I next discuss the violation of yet 
another classical inequality namely, the Cauchy-Schwarz inequality. 

Violation of the Cauchy-Schwarz Inequality 

The other signature for the non-classical nature of the emitted fields is the 
violation of Cauchy-Schwarz inequality [6] . The degree of violation is determined 
by the ratio R = (9^31 (t))'^ / g^^ {0)g^i (t), i? < 1 for a classical source. One 
needs to optimize the emissions in such a way that one achieves a large i?™"^. 
In the present scheme, this condition essentially boils down to demanding sub- 
poissonian statistics in the Qa {t) for ? = 1,3 -which ensures that the photons 
are emitted at spaced out intervals - in addition to large values for ^g-^ (t). The 
function gl^ (r) can be tailored to show sub-poissonian statistics by making 
the coupling field weak { Vli < r2) . The auto-correlation function 333 (t) 
shows sub-poissonian statistics even for ^3 > f 4 since the population of the 
excited levels is smaller. Fig4 shows the correlation functions and the ratio R 
for various parameter values. Here, the dimensionless variable t = ^t where 

(2) 

7 = 2itMHz. The function .931 (t) increases with fi^/ (Fig4 (c)), while both 

the auto-correlation functions decrease with fi^/ (Fig4 (a) and (b)). On the 

f 2) 

other hand, the time delay for the occurrence of peak value of g\^ {t) is smaller 
for larger values of il^/. This means that, in principle, the occurrence of i?™"^ 
can be delayed/advanced by decreasing/increasing the coupling field strengths 
- implying a trade-off between the value of R™-"-^ and the time delay for the 
occurrence of the same. Note that R"^"-^ = IQ^ — 10^ as shown in Fig4 (d) and 
increases with $7^/ . It is noteworthy that the oscillatory behaviour of R for 
large Vlrf provides different time delays at which the violation is maximum. 
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Lastly, the merit in this scheme is that the the absorption profile being very 
narrow (of the order of the natural hne- width) [19j . the excitation/deexcitation 
of the g.s. and the upper excited state occurs within a very narrow frequency 
band. This is because of the three photon resonance which induces atomic 
coherence. It has been predicted earher [11] that when a three-level cascade 
system showing EIT is coupled to a third driving field, the ground state (g.s.) 
shows a very narrow absorption within the EIT window at three photon reso- 
nance; the Doppler integrated absorption does not suffer much change in this 
geometry. This type of narrow absorption features and three-peak absorption, 
in four- level systems, have been subsequently reported experimentally [12]. An 
insight into this behavior is revealed by the study of the single atom dynamics 
[2U] : when the middle transition coupling is strong {^rf in this case) the four 
level atom behaves like a two level atom spanned by levels |1 > and |4 > hence 
inducing a strong atomic coherence between these two levels. A steady state 
analysis further reveals that in an ensemble this feature manifests as a dom- 
inant contribution of the third order nonlinear susceptibility (due to the strong 
atomic coherence) as compared to the linear susceptibility. This third order 
nonlinear susceptibility is responsible for the efficient transfer of population to 
the uppermost level by inducing a narrow absorption in the ground state. The 
role of rirf is not only to transfer the population to level |4 > but the strong 
atomic coherence it introduces is responsible for the strong correlation between 
the photons emitted by the extreme transitions. 

The absorption features for the model under consideration is very similar, 
and since the absorption occurs within an EIT window, there is an absorption 
free zone [TH] for detunings larger than the line width of the central peak, and 
detunings less than the FWHM of the EIT window, i.e. for 15MHz > |Ai| > 
6MHz. This ensures cleaner signals and a reduction in the noise due to the 
absorption and emission into other modes close to the resonant frequency which 
is a major advantage in this system. I would like to mention here that if the 
coupling between the levels |2 > — 13 > were in the visible region instead of radio- 
frequency region then the absorption of the |3 > — 14 > transition diminishes due 
to Doppler averaging [11]. Even though this does not affect the antibunching 

(2) 

nature of g^i (t), the strength of the correlation function would become weaker. 
Hence the signature for the violation of the Cauchy-Schwarz inequality may also 
become weaker. 

To summarize, I propose a model for generating non-classical pairs, wherein 
Doppler effects are negligible at room temperatures, and also the absorption/emission 
proceeds in a narrow frequency band . The violation of the Cauchy-Schwarz in- 
equality is orders of magnitude larger than the experimentally reported data, so 
far, and therefore could provide very clear signals in the near future. Another 
interesting feature is the antibunching nature of the cross-correlation between 
the emission of photons from the extreme transition which occurs because of 
the cancellation of emission due to different two-photon pathways. Thus the 
photon pair generated in this scheme shows the violation of two of the classical 
inequalities. 
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Appendix 

The solutions for the density matrix elements in the Laplace space for various 
initial conditions and both the strong coupling of flrf and weak fl^f are listed 
below. 

Strong flj-f 

For the initial condition p(0) = |3 >< 3|, the relevant Pij{T) in the Laplace 
space is given by 

C.{P22{t)) = i)7{s) 

= + f 2 + f 3 + 2ni)(l + 5^(1 - S 

-fi){s + f 2 + f4)) + ^0202(1 _ s _ fa)] 

where = Uf^iis + f 4 + aj), ^3 = nf^3(s + a^), ^2 = ^t=i{s + ttj) 
>C(p33(r)) = Ms) 

= i[(s'+fi+f2f3+,s(f3+2f2) 

+2f^i)(i + ;iiS;)(- + r2 + r4) 

For the initial condition p(0) = |2 >< 2| we have 



^ + f 2 + s(f 2 + 2f 3) + f 2r3 + 



2iil) + ^^{2nlnUs - 1) - 2r!2(s + f3) 

(s2 + f i + s(f 2 + 2f 3) + 2nl + t^fll))] 
For the initial condition p(0) = |1 >< 1| we get 
C.{P22{t)) = ^s) 

= »["i + (« + r3)((s + f3 + f2) 

(s + f3) + 0i)] 

Weak Vlrf 

For the initial condition p(0) = |3 >< 3| wc have 

C{p22{t)) = -07(5) 
1 anl 



= ik;^s+ ^^Mi - 2(s + T2)n2C2 + 
+ Q2^3{{S + r3)C3 - 2O3C2)) 



15 



-{Tl + s{s - 20f) + 2r2(s - nl)) 



((1 + 202C2)(s2 + f 1 + f 3(S + f4) + 

2r4S + 2Qi) + 2C302r23(s + f 4 - 1))] 



-^[(1 + 2r!2C2)Gs2 + fl + f3(s + f4) 



2^1 + 2sT4,) + 2n2^3C3(s + 2r4 + 1)] 



= + 2O2C2) - 02C3(s + r4)] 

For the initial condition p(0) = |2 >< 2| we have 

jO.{P22{t)) = -07(5) 

= 3fe[?+2^^i"2Ci-(. + f2)(l- 

2fl^C2) + ^^^^{{{s + f 3)C3 - 2O3C2)) 
+ (5Tr^(ri + .(«-20?)+2f2(«-f2?) 



(2f]2C2(s2 + f 1 + f 3(5 + f 4) + 2r4S 
+202) ^ 2C3fi2^^3(s + f 4 - 1)))] 

For the initial condition p(0) = |1 >< 1| 

£(p22(r))= V7(s) = 20f/(srf2p) 



where 



Ci= -f^if72[(s + f4)(s + f4 + 

f2) + 0?-17i]/d4p 
C2= -02[(s + f4)(s + f3)(s + f4 + f2) + 
fllis + f 2 + f 4) + (S + f 3)0?)]/rf4p 
C3= 0203[-(s + f3)(s + f4) + 

Here (i4p = nf^g(s + a,), ^3^ = nf^3(s + Ui) and d2p = (s + f 2)^ + 4f2?. 
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